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Abstract. We describe odd isotopes as formed of a core plus one neutron. We have calculated the modifica- 
tion of single neutron energies brought by couplings of the neutron with collective 2"'" phonons in the cores. 
The results reproduce very well the inversion of 2s and 1^5/2 shells in carbon isotopes up to ^^C while in 
oxygen isotopes the correction is also large but do not show any inversion in agreement with experiments. 
The calculated energies are close to the experimental ones in both series of isotopes except in ^^C for the 
2s state for which our coupling is too weak. 
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1 Introduction 



In the last decade the properties of carbon and oxygen iso- 
topes have been extensively studied experimentally and 
theoretically. They show very interesting new behaviour 
when one compares with heavier nuclei. In particular car- 
bon isotopes offer a very interesting situation. Indeed in 
lighter carbons, including ^^C, one notices an inversion of 
the 2s and 1^5/2 shells, inversion which disappears when 
N, the number of neutrons, is equal to 12 or 14. At con- 
trary, in oxygen isotopes the normal sequence of 1^5/2 and 
2s shells with the 1^5/2 lower than the 2s one is observed 
but neutron energies for states closed to the Fermi surface 
are far from energies given by a Saxon Woods potential or 
an Hartree Fock potential P] . Therefore we are interested 
to see what is the origin of this inversion in carbon nuclei 
and how the situation is different in oxygen nuclei. 

These nuclei have been the object of calculations mostly 
using shell models but we present a different description 
based on the assumption of a two-body system formed 
of a neutron added to an inert core for which the two- 
body correlations are introduced via a coupling between 
the neutron and a phonon of the core. The corrections to 
the one-body neutron-core potential [^ due to this type of 
coupling induces a modification of the single neutron ener- 
gies which can be easily calculated from the known prop- 
erties of the phonons of the core. In section [5] we present 
our model and detail the calculation of corrections due 
to two-body correlations introduced via neutron-phonon 
couplings. Then in section [3] we present and discuss the 
corresponding results obtained for carbon and oxygen iso- 
topes. At last section m is devoted to our conclusion. 



2 Description of isotopes with an odd 
number of neutrons and contribution of 
Neutron-phonon couplings to the single 
neutron energies 

Assuming a neutron plus an inert core, the neutron states 
will be calculated first using a Saxon- Woods one-body po- 
tential for the neutron-core interaction written as [5]: 



V{r) 



Vq fir) - 0.44r2(l.s) 



where 



f{r) = 



1 + exp 



r-Rr 



1 dfjr) 
r dr 



and 



(1) 



(2) 



Vb = (50.5-32.5 ^ ^ j MeV (3) 

with a parameter a of 0.75 fm, larger than usual, to take 
account of the diffuse surface of such light nuclei. The 
other parameters are Rq = roA^'^ with ro=1.27 fm where 
A is the nucleon number in the core. These parameters 
were fitted in order to reproduce at best the neutron shells 
in nuclei with ^ > 40 where particle-phonon couplings are 
expected to have little effect on single particle energies, 
therefore it can be considered as a good phenomenological 
representation of the first order one-body potential. 

To this first order potential we add the contribution 
due to two-body correlations taken into account via a cou- 
pling between the neutron and a phonon of the core. The 
corresponding modification of the potential is given in ref- 
erences minis]- To calculate the corresponding modifica- 
tion to the single neutron energies we follow reference [5] . 
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This correction, Z\e„ for the neutron in state n is given 
by: 



Z\e„ = En_xFj 



N,X^N,\ 



(t>n{'r)VoN{r)(l)\{r)dr 



(4) 



where: 

- N and A mean respectively the phonons which are in- 
cluded in the calculation and the state of the intermediate 
neutron. 

- (/'n(r) and 4'\{v) are respectively the wave functions 
calculated in the potential of Eq. [T] for the neutron in 
states n and A with respective energies e„ and eA- 

- VoAT is the amplitude for exciting the phonon N of 
angular momentum L. As in reference [2] we make a phe- 
nomenological approach and write Vqat as: 



Von{v) 



1 



(2L + 1)1/2 



PlR^ 



dU{r ). 
dr 



r-*K) (5) 



where U(r) is the central part of our Saxon Woods poten- 
tial of Eqs. (fT]l3)) and Pl is related to the B(EL) by: 






B{EL) 
Fat A is a coefficient given by: 
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(7) 



where n\ is the occupation number of the state A and En 
is the energy of the studied neutron after correction, which 
is then given by: 



En — £« + ^£n 



(8) 



These equations have been solved by iteration to get 
the new energy E„. 



3 Results for carbon and oxygen isotopes 

3.1 Carbon isotopes 

All the results are given in table [T] for the states closed 
to the Fermi surface. We first concentrate on the ener- 
gies of the 2s and 1^5/2 shells which present an unusual 
behaviour. We first calculate their energies in our Saxon 
Woods potential of Eq. [ij We see that up to ^°C the 2s 
state is lower than the 1^5/2, then the order of the two 
shells is inverted and the usual situation is recovered for 
^''C. The same situation was already observed for ^^Be 
in [5] when a pure Saxon Woods potential was used. This 
property is clearly related to the radius of the neutron-core 
potential. Indeed for smaller radii the centrifugal barrier is 
more efficient and pushes up the d-states compared to the 
2s one. However the calculated energies are quite far from 
the experimental ones pr7ll51[giPIIl[TT|[ni[ni[TilP^[TCl[r7l 
[T51I19U201I2I1I221I23) which are also given in the table H 



Table 1. Results for carbon isotopes and comparison with 
experimental values when they are known. 







2Si/2 


ld5/2 


1^3/2 


2d5/2 


2rf3/2 




En 


-1.02 


-0.76 


1.80 


1.94 


3.83 


12^ 


A€n 


-0.93 


-0.38 


-0.04 


-0.08 


-0.67 


E„ 


-1.95 


-1.14 


1.76 


1.86 


3.16 




-LJexp 


-1.86 


-1.12 




1.74 


3.3 




fn 


-0.93 


-0.73 


1.76 


1.93 


3.6 


"C 


A€n 


-0.21 


-0.06 


-0.01 


-0.01 


-0.16 


E„ 


-1.13 


-0.79 


1.75 


1.92 


3.44 




-LJexp 


-1.22 


-0.48 






3.47 




En 


-0.90 


-0.78 


1.72 


1.92 


3.38 


16^ 


Ae,, 


-0.37 


0. 


-0.03 


0.03 


-0.11 


E„ 


-1.27 


-0.78 


1.69 


1.95 


3.27 




t^exp 


-0.73 


-0.48 






3.2 




Cn 


-0.93 


-1.02 


1.63 


1.89 


2.91 


20 c 


A€„ 


0.37 


-0.18 


0.05 


0.01 


-0.02 


E„ 


-0.56 


-1.20 


1.68 


1.90 


2.89 




Eexp 


0.5 








3. 



Table 2. Energies and B{E2) of the 2^ 
calculation for the carbon cores 



states used in our 





12c 


"C 


16(. 


20 c 


E2 (MeV) 


4.43 


7.01 


1.76 


1.62 


B(E2) (e^fm^) 


39.7 


18.7 


13. 


18.4 



Therefore we expect two body correlations to be impor- 
tant as it was the case for the inversion of lpi/2 and 2s 
shells in i^Be, ^^Li and ^'^Be. Indeed the cores that we are 
considering have 2"*" states with a low excitation energy 
and very large B{E2) due to a large collectivity. Such col- 
lective 2+ states were shown to be responsible for a large 
modification of neutron energies in a core of ^^Be and to 
give most of the contribution to the corrective term [51 
[21] • We give in Table II the energies and B(E2) for the 
2"*" excited states that we have included in our calcula- 
tion [ 25] , assuming that it gives most of the contribution. 
The results for cores of ^^C, ^^C, ^^C and ^"C are given 
in the table [T] where we see that now the inversion of the 
two shells is amplified and individual neutron energies for 
2s and 1^5/2 shells are close to the experimental ones. For 
a i^C-core we did not calculate the corrective term. In- 
deed from the calculation in a SW potential and as given 
by experimental energies the two shells are very closed, so 
that i^C is very likely a mixture of (2s)^ and (1^5/2)^ con- 
figurations for the two last neutrons. It is indeed verified 
when it is described as a core of ^^C plus two neutrons [26j . 
Therefore in our model it may not be taken as a core. 

The same calculation has been done for 1^3/2, 2^5/2 
and 2^3/2 shells and results are given in table [T] For higher 
neutron states the corrections due to the couplings are 
negligible as expected. Our 2^5/2 reproduces very well the 
second 5/2+ excited state in ^^C [T5^ while there is no 
experimental 3/2"*" state corresponding to our 1^3/2 state. 
This state is very close to the second 5/2+ and perhaps 
not easily detected experimentally. However our 2^3/2 is 
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Table 3. Results for oxygen isotopes and comparison with 
experimental values when they are known [28II30IIM] . 



Table 4. Energies and B (£2) of the 2^ 
cores used in our calculation. 



states in the oxygen 







2Sl/2 


ld5/2 


1^3/2 




En 


-3.5 


-4.92 


1.53 


14q 


A€n 


-0.89 


-0.79 


-0.27 


E„ 


-3.90 


-5.02 


0.99 




Eexp 


-4.50 


-6.4 


1.54 




f-n 


-3.2 


-4.52 


1.37 


16q 


Af..a 


-0.61 


-0.54 


-0.5 


E„ 


-3.81 


-4.05 


0.87 




Eexp 


-3.3 


-4.14 


0.94 




Cn 


-2.85 


-4.1 


0.65 


220 


Atn 


0.12 


0.15 


0.26 


E„ 


-2.73 


-3.95 


0.91 




Eexp 


-2.74 


-3.8 






En 


-2.82 


-4.08 


0.38 


24q 


Aen 
E„ 










Eea;p 


-4.09 


-5.05 


0.77 



in great agreement with the known high 3/2^ state in all 
our systems even though the corrective term is very small. 

For a ^°C core we recover the usual order with the 
1(^5/2 lower than the 2s, therefore ^°C may be assumed 
to have a closed 1^5/2 shell. However for the 2s state the 
neutron-phonon coupling is not strong enough to make the 
state unbound as it should. In this case we may expect a 
further contribution to the potential of a term correspond- 
ing to the excitation of collective pairing vibrations in the 
N-2 and iV-|-2 systems g], therefore in ^^C and '^'^C. In- 
deed in particular in ^^C it has been shown that the wave 
functions are mixture of several configurations for the two 
last neutrons |26j . This type of contribution was found 
long ago to have a noticeable contribution in ''"Ca [27] . 

As a summary we can say that, even though we have 
restricted the neutron-phonon coupling to the low energy 
2+ phonon, we have been able to reproduce the inversion 
of 2s-l(i5/2 shells in the lightest carbon isotopes (up to 
^"^C) and could get reasonable agreement with experimen- 
tal neutron energies for all calculated states. 



3.2 Oxygen isotopes 

We have performed the same calculations for oxygen iso- 
topes assumed to be described as an inert core plus one 
neutron. We have considered cores of ^''O, ^^O, ^°0 and 
^^O.^^O may not be assumed as a core. Indeed because of 
the proximity of the 2s and 1^5/2 shells, the wave func- 
tion is a mixture of two configurations with the two last 
neutrons in one of these two shells. As for carbons we 
have calculated the contribution to the neutron energies 
of neutron-phonon couplings due to the low energy 2+ ex- 
cited state of the core. The energies and B(E2) are taken 
from measurements 28 , 29 and given in table |H The re- 
sults of our calculation are reported in table[3]for 2s, 1^5/2 
and 1^3/2 neutron states together with the experimental 
energies [SHllSnilSIllSlES] ■ 





14q 


16(3 


220 


E2 (MeV) 


6.59 


6.92 


3.19 


B(E2) (e^fm^) 


43. 


40.6 


21. 



We see first that we have no inversion of the 2 s and 
1^5/2 shells for any of the isotopes in agreement with 
experiments. On the other side the neutron-phonon cou- 
plings are quite important and put the calculated energies 
close to the experimental ones. As in carbon isotopes, for 
higher states and already for the 2^5/2 and 2^3/2 states, 
these corrections are negligible (<0.1 MeV). For ^'^O we 
do not know any 2+ state then could not calculate the 
coupling term. However we note that to obtain the exper- 
imental energy, in particular for the 2s state, we need a 
large corrective term what suggests that there is a low 2+ 
state in ^^O with a large B(E2). This is plausible since in 
^■^O the l(i5/2 and 2s shells are closed to each other and 
a 2+ state can be constructed by coupling two neutrons 
of the 1^5/2 shell to J'^ = 2+. The excitation energy of 
this state should be then close to the low 2+ state of ^^O 
known at 1.98 MeV. 



4 Conclusions 

To conclude we were able to understand and reproduce 
the properties of single neutron states in carbon and oxy- 
gen isotopes with an odd number of neutrons. We have 
shown the importance of two body correlations which are 
treated in our work as neutron-phonon couplings reduced 
to coupling with the lowest 2+ state in the core. In par- 
ticular we could explain that for carbon isotopes the 2s 
and l(i5/2 shells are inverted up to ^"C as shown exper- 
imentally. For oxygen isotopes no such inversion appears 
in our calculations as required by experiments. For both 
series of nuclei we get neutron energies closed to exper- 
imental values except for ^^C where the coupling of the 
neutron with the 2+ phonon of ^"C is not strong enough 
to get the measured unbound energy of the 2s state even 
though it reduces the discrepancy between calculation and 
measurement. 

Apart the corrections on single neutron energies, particle- 
phonon couplings induce a modification of the wave func- 
tions which are now a mixture of one single state with con- 
figurations where one neutron is coupled to phonons [2]. 
One, sometimes, makes a parallel between this mixture 
of complex configurations and a deformation of the cor- 
responding state as it is explicitly used in sulfur stud- 
ies [3S]. This equivalence was already clear in the case of 
-'^^Be where the inversion of 1/2+ and l/2~ was shown to 
be due to two-body correlations in a calculation similar 
to the present one E] or to a deformation of the one-body 
average potential j36ll37j due to the presence of the low 2+ 
state in ^°Be which was considered as a vibrational state 
in the first case and a rotational state in the second one. 
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